Hadronic Matter is Soft 
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The stiffness of the hadronic equation of state has been extracted from the production rate of 
K+ mesons in heavy ion collisions around 1 A GeV incident energy. The data are best described 
with a compressibility coefficient k around 200 MeV, a value which is usually called "soft". This 
is concluded from a detailed comparison of the results of transport theories with the experimental 
data using two different procedures: (i) the energy dependence of the ratio of K + from Au+Au and 
C+C collisions and (ii) the centrality dependence of the K + multiplicities. It is demonstrated that 
input quantities of these transport theories which are not precisely known, like the kaon-nucleon 
potential, the AN — > NK + A cross section or the life time of the A in matter do not modify this 
conclusion. 
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Since many years one of the most important challenges 
in nuclear physics is to determine E / A(p,T), the en- 
crgy/nucleon in nuclear matter in thermal equilibrium 
as a function of the density p and the temperature T. 
Only at equilibrium density, po, the energy per nucleon 
E/A(p — po,T = 0) = —16 MeV is known by extrapo- 
lation of the Weizsacker mass formula to infinite matter. 
This quest has been dubbed "search for the nuclear equa- 
tion of state (EoS)" . 

Modelling of neutron stars or supernovae have not yet 
constrained the nuclear equation of state Therefore, 
the most promising approach to extract E/A(p,T) are 
heavy ion reactions in which the density of the colliding 
nuclei changes significantly. Three principal experimen- 
tal observables have been suggested in the course of this 
quest which carry - according to theoretical calculations 
- information on the nuclear EoS: (i) the strength distri- 
bution of giant isoscalar monopole resonances @, 0] , (ii) 
the in-plane sidewards flow of nucleons in semi-central 
heavy ion reaction at energies between 100 A MeV and 
400 A MeV 0] and (iii) the production of K+ mesons m 
heavy ion reactions at energies around 1 A GeV Al- 
though theory has predicted these effects qualitatively, a 
quantitative approach is confronted with two challenges: 
a) The nucleus is finite and surface effects are not neg- 
ligible, even for the largest nuclei and b) in heavy ion 
reactions the reacting system does not come into equilib- 
rium. Therefore complicated non-equilibrium transport 
theories have to be employed and the conclusion on the 
nuclear equation of state can only be indirect. 

(i) The study of monopole vibrations has been very 
successful, but the variation in density is minute. There- 
fore, giant monopole resonances are sensitive to the en- 
ergy which is necessary to change the density of a cold nu- 
cleus close to the equilibrium point pq. According to the- 
ory the vibration frequency depends directly on the force 
which counteracts to any deviation from the equilibrium 
and therefore to the potential energy. The careful analy- 
sis of the isoscalar monopole strength in non-rclativistic 



and relativistic mean field models has recently con- 
verged Q due to a new parametrization of the relativis- 
tic potential. These calculations allow now for the deter- 
mination of the compressibility n = 9p 2 d E ^ p p ' T ^ \p=p 
which measures the curvature of E/A(p,T) at the equi- 
librium point. The values found are around k = 240 
MeV and therefore close to what has been dubbed "soft 
equation of state" . 

(ii) If the overlap zone of projectile and target becomes 
considerably compressed in semi-central heavy-ion colli- 
sions, an in-plane flow is created due to the transverse 
pressure on the baryons outside of the interaction region 
with this flow being proportional to the transverse pres- 
sure. In order to obtain a noticeable compression, the 
beam energy has to be large as compared to the Fermi 
energy of the nucleons inside the nuclei and hence a beam 
energy of at least 100 A MeV is necessary. Compression 
goes along with excitation and therefore the comprcs- 
sional energy of excited nuclear matter is encoded in the 
in-plane flow. It has recently been demonstrated [|| that 
transport theories do not agree quantitatively yet and 
therefore former conclusions j?} have to be considered as 
premature. 

(iii) The third method is most promising for the study 
of nuclear matter at high densities and is subject of this 
Letter. K + mesons produced far below the NN thresh- 
old cannot be created in first-chance collisions between 
projectile and target nucleons. They do not provide suf- 
ficient energy even if one includes the Fermi motion. The 
effective energy for the production of a K + meson in the 
NN center of mass system is 671 MeV as in addition 
to the mass of the kaon a nucleon has to be converted 
into a A to conserve strangeness. Before nucleons can 
create a K + at these subthreshold energies, they have to 
accumulate energy. The most effective way to do this is 
the conversion of a nucleon into a A and to produce in 
a subsequent collision a K + meson via AiV — * NK + A. 
Two effects link the yield of produced K + with the den- 
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sity reached in the collision and the stiffness of the EoS. 
If less energy is needed to compress matter (i) more 
energy is available for the K + production and (ii) the 
density which can be reached in these reactions will be 
higher. Higher density means a smaller mean free path 
and therefore the A will interact more often increasing 
the probability to produce a K + and hence, it has a 
lower chance to decay before it interacts. Consequently 
the K + yield depends on the comprcssional energy. At 
beam energies around 1 A GeV matter becomes highly 
excited and mesons are formed. Therefore this process 
tests highly excited hadronic matter. At beam energies 
> 2 A GeV first-chance collisions dominate and this sen- 
sitivity is lost. 

In this Letter we would like to report that for the third 
approach different transport theories have converged. 
Two independent experimental observables, the ratio of 
the excitation functions of the K + production for Au+Au 
and for C+C and a new observable, the depen- 

dence on the number of participants of the K + yield show 
that nucleons interact with a potential which corresponds 
to a compressibility of n < 200 McV in infinite matter in 
thermal equilibrium. This value extracted for hadronic 
matter at densities around 2.5 times the normal nuclear 
matter density is very similar to that extracted at normal 
nuclear matter density. A key point of this paper is to 
demonstrate that the different implementation of yet un- 
solved physical questions, like the NA — > K + AN cross 
section, the KN interaction as well as the life time of the 
nuclear resonances in the hadronic environment do not 
affect this conclusion. 

In order to determine the energy which is necessary 
to compress infinite nuclear matter in thermal equilib- 
rium by heavy ion reactions in which no equilibrium is 
obtained one chooses the following strategy: The trans- 
port theory calculates the time evolution of the quantal 
particles described by Gaussian wave functions. The time 
evolution is given by a variational principle and the equa- 
tions one obtains for this choice of the wave function are 
identical to the classical Hamilton equations where the 
classical two-body potential is replaced by the expecta- 
tion value of the real part of the Bruckner G-matrix. For 
this potential the potential energy in infinite nuclear mat- 
ter is calculated. To determine the nuclear equation of 
state we average this (momentum-dependent) two-body 
potential over the momentum distribution of a given tem- 
perature T and add to it the kinetic energy. Expressed 
as a function of the density we obtain the desired nuclear 
equation of state E/A(p, T). Our two-body potential has 
five parameters which are fixed by the binding energy of 
infinite nuclear matter at po, the compressibility k and 
the optical potential which has been measured in pA re- 
actions. 

Once the parameters are fixed we use the two-body 
potential with these parameters in the transport calcula- 
tion. There is an infinite number of two-body potentials 
which give the same equation of state because the range 
of the potential does not play a role in infinite matter. 



The nuclear surface measured in electron scattering on 
nuclei fixes the range, however, quite well. The uncer- 
tainty which remains is of little relevance here (in con- 
tradiction to the calculation of the in-plane flow which is 
very sensitive to the exact surface properties of the nuclei 
and hence to the range of the potential). 

We employ the Isospin Quantum Molecular Dynam- 
ics (IQMD) Q approach with the following equations of 
motion: 



and f.-m. 

Pi ~ dn and n ~ dp, 



(1) 



where the expectation value of the total Hamiltonian 
reads as (H) = (T) + {V) with 

i T ) = yiL 

X 1 ^ 2m, 

< y > = T,/Zj 'fi(r,P,t)V^ <fj{r ',p',t)drdr'dpdp\2) 

i j>i 

and fi being the Gaussian Wigner density of nucleon 
i. The baryon-potential consists of the real part of 
the G-Matrix which is supplemented by the Coulomb 
interaction between the charged particles. The former 
can be further subdivided in a part containing the con- 
tact Skyrme-type interaction only, a contribution due 
to a finite range Yukawa-potential, and a momentum- 
dependent part with 
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with Zi, Zj the charges of the baryons i and j. For more 
details we refer to Ref. Q- 

We include in this calculation all inelastic cross sections 
which are relevant for the K + production. For details 
of these cross sections we refer to ^} . Unless specified 
differently, the change of the K + mass due to the kaon- 
nucleon (KN) interaction according to m K (p) = mff(l — 
0.075—) is taken into account, in agreement with recent 
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self-consistent calculations of the spectral function of the 
K + ^jj- The A potential is 2/3 of the nucleon potential, 
assuming that the s quark is inert. 

In order to minimize the experimental systematical er- 
rors and the consequences of theoretical uncertainties it is 
better to compare ratios of cross sections rather than the 
absolute values 01 ■ We have made sure that the stan- 
dard version of IQMD reproduces the excitation function 
for Au+Au as well as for C+C quite well [l^. These ra- 
tios are quite sensitive to the nuclear potentials because 
the compression obtained in the Au+Au collisions is con- 
siderable (up to 3po) and depends on the nuclear equa- 
tion of state whereas in C+C collisions the compression 
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is small and almost independent on the stiffness of the 
EoS. 

Figure ^ shows the comparison of the measured ratio 
of the K + multiplicities obtained in Au+Au and C+C 
reactions together with transport model calculations 
as a function of the beam energy. 

We see clearly that the form of the yield ratio depends 
on the potential parameters (hard EoS: k = 380 MeV, 
thin lines and solid symbols, soft EoS: k — 200 MeV , 
thick lines and open symbols) in a quite sensible way 
and that the prediction in the standard version of the 
simulation (squares) for a soft and a hard EoS potential 
differ much more than the experimental uncertainties. 
The calculation of Fuchs et al. [3| given in the same 
graph, agrees well with our findings. 

This observation is, however, not sufficient to deter- 
mine the potential parameters uniquely because in these 
transport theories several not precisely known processes 
are encoded. Therefore, it is necessary to verify that 
these uncertainties do not render this conclusion prema- 
ture. Figure ^ top, shows as well the influence of the 
unknown NA — > K + AN cross section on this ratio. We 
confront the standard IQMD option (with cross sections 
for AN interactions from Tsushima et al. ^(J) with an- 
other option, er(AA) = 3/4a(NN) which is based 
on isospin arguments and has been frequently employed. 
Both cross sections differ by up to a factor of ten and 
change significantly the absolute yield of K + in heavy 
ion reactions but do not change the shape of the ratio. 

The middle part demonstrates the influence of the 
kaon-nucleon potential which is not precisely known at 
the densities obtained in this reaction. The uncertain- 
ties due to the A life time are discussed in the bottom 
part. Both calculations represent the two extreme val- 
ues for this lifetime [Tol ] which is important because the 
disintegration of the A resonance competes with the K + 
production. 

Thus we see that these uncertainties do not influence 
the conclusion that the excitation function of the ratio is 
quite different for a soft EoS potential as compared to a 
hard one and that the data of the KaoS collaboration are 
only compatible with the soft EoS potential. The only 
possibility to change this conclusions is the assumption 
that the cross sections are explicitly density dependent 
in a way that the increasing density is compensated by a 
decreasing cross section. It would have a strong influence 
on other observables which are presently well predicted 
by the IQMD calculations. 

We would like to add that the smoothness of the excita- 
tion function also demonstrates that there are no density 
isomers in the density regions which are obtained in these 
reactions because the K + excitation function would be 
very sensitive to such an isomeric state (l6j . 

The conclusion that nuclear matter is best described by 
a soft EoS, is supported by another variable, the depen- 
dence of the K + yield on the number of participating nu- 
clcons Ap ar t ■ The prediction of the IQMD simulations in 
the standard version for this observable is shown in Fig. [21 
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FIG. 1: Comparison of the measured excitation function of 
the ratio of the K + multiplicities per mass number A obtained 
in Au+Au and in C+C reactions (Ref. |T^ with various cal- 
culations. The use of a hard EoS is denoted by thin (blue) 
lines, a soft EoS by thick (red) lines. The calculated energies 
are given by the symbols, the lines are drawn to guide the eye. 
On top, two different versions of the NA — > K + AN cross sec- 
tions are used. One is based on isospin arguments IHi . the 
other is determined by a relativistic tree level calculation lrH| . 
The calculation by Fuchs fl^l are shown as dotted lines. Mid- 
dle: IQMD calculations with and without KN potential are 
compared. Bottom: The influence of different options for the 
life time of A in matter is demonstrated. 



The top of the figure shows the kaon yield M K + /A part 
for Au+Au collisions at 1.5 A GeV as a function of the 
participant number A part for a soft EoS using different 
options: standard version (soft, KN), calculations with- 
out kaon-nucleon interaction (soft, no KN) and with the 
isospin based NA — > AAA' + cross section (soft, KN, 
a*). These calculations arc confronted with a standard 
calculation using the hard EoS potential. The scaling of 
the kaon yield with the participant number can be pa- 
rameterized by M K + = Ap art . 

All calculations with a soft EoS show a rather simi- 
lar value of a - although the yields are very different - 
while the calculation using a hard equation shows a much 
smaller value. Therefore we can conclude that also the 
slope value a is a rather robust observable. 

The bottom of Fig. |3 shows that a depends smoothly 
on the compressibility k of the EoS. Whether we in- 
clude the momentum dependence of the nucleon nucleon 
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FIG. 2: Dependence of the K + scaling on the nuclear equation 
of state. We present this dependence in form of M K + = Ap art . 
On the top the dependence of M K + /Ap art as a function of 
A part is shown for different options: a "hard" EoS with KN 
potential (solid line), the other three lines show a "soft" EoS, 
without KN potential and a(NA) from Tsushima fl^l (dot- 
ted line), with KN potential and the same parametrization 
of the cross section (dashed line) and with KN potential and 
a(NA) = 3/4a(NN). On the bottom the fit exponent a 
is shown as a function of the compressibility for calculations 
with momentum-dependent interactions (mdi) and for static 
interactions (t^ = 0, dashed line). 



interaction (with mdi) or not (without mdi) does not 
change the value of a as long as the compressibility is not 
changed - in stark contrast to the in-plane flow. Again, 
the measured centrality dependence for Au+Au at 1.5 A 
GeV from the KaoS collaboration [j^. a = 1.34 ± 0.16, 
is only compatible with a soft EoS potential. 

This finding is also supported by a more recent anal- 
ysis 0, |2(| of the in-plane flow which supersedes the 
former conclusion that the EoS is hard |2l| (made before 
the momentum-dependent interaction has been included 
in the calculations). Due to the strong dependence of 
the in-plane flow on the potential range parameter and 
its dependence on the particles observed these conclu- 
sions are much less firm presently. Comparisons of the 
out-of-plane squeeze of baryons also show a preference 
for a soft equation of state with momentum dependent 
interactions (23. 

In conclusion, we have shown that the two experimen- 
tal observables which are most sensitive to the potential 
parameters of the nucleon-nucleon interaction are only 
compatible with those parameters which lead in nuclear 
matter to a soft hadronic EoS. This conclusion is robust. 
Uncertainties of the input in these calculations, like the 
KN potential at high densities, the lifetime of the A in 
matter and the A7V — > NK + A cross section do not influ- 
ence this conclusion. The potential parameter k is even 
smaller than that extracted from the giant monopole vi- 
brations. Thus the energy which is needed to compress 
hadronic matter of k < 200 MeV is close to the lower 
bound of the interval which has been discussed in the 
past. 

We would like to thank all members of the KaoS Col- 
laboration for fruitful discussions especially A. Forster, 
P. Senger, C. Sturm, and F. Uhlig. 
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